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We have measured the differential cross section for the inclusive production of ip(2S) mesons 
decaying to [i + fi~ that were produced in prompt or B-decay processes from pp collisions at 1.96 
TeV. These measurements have been made using a data set from an integrated luminosity of 1.1 fb _1 
collected by the CDF II detector at Fermilab. For events with transverse momentum pt(V'(2S)) > 
2 GeV/c and rapidity |y(i/)(2S))| < 0.6 we measure the integrated inclusive cross section o~(pp — > 
tp{2S)X) ■ Br(^(2S) -> to be 3.29 ± 0.04(stat.) ± 0.32{syst.) nb. 



The mechanism for producing heavy vector mesons in pp collisions is not well understood. The experimental 
measurement of prompt 3/ip and ip(2S) production cross sections by CDF in Tevatron Run I [l| showed that the 
measured cross sections were one to two orders of magnitude larger than expected from the leading order (LO) 
color-singlet models. Theoretical efforts to improve the calculations added color octet contributions that increased 
the predicted cross sections, e.g., in the non-relativistic QCD model (NRQCD) [2]. Recently there have been other 
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approaches that do not directly introduce a color octet amplitude; rather, they incorporate the effects of multiple 
gluon processes during the production process, e.g., the fc^-factorization formalism Q and the gluon tower model Q|. 

The NRQCD model with parametrized production matrix elements adjusted to data can successfully account for the 
Tevatron prompt ^(2S) cross section measurements, but it makes an unequivocal prediction of increasing transverse 
polarization of vector mesons as their transverse momentum px from production increases Q • A recent polarization 
measurement at CDF Q contradicts the NRQCD model prediction. 

Experimentally, the extraction of direct J ftp production information is complicated by significant feed-down from 
decays of promptly-produced higher-mass charmonium states (xo V*(2S)) to 3/ip mesons. This is not a problem for 
direct ip(2S) production because there are no reported charmonium states with significant hadronic production cross 
sections that decay to the ip(2S). Consequently the ip(2S) provides an ideal testing ground for studying charmonium 
hadroproduction mechanisms. In this paper we present a measurement of the pr dependence of the V>(2S) production 
cross section over the ip(2S) transverse momentum range 2 < pr(-0(2S)) < 30 GeV/c with rapidity | y(ip(2S)) |< 0.6. 
This measurement greatly increases the statistical power of the data in the perturbative regime (mr 3> Aqcd), 
facilitating comparison with theory. 

We use data taken using the CDF II detector at the Fermilab Tevatron at 1.96 TeV 0]. The integrated luminosity 
of the data sample is 1.1 fb -1 . The CDF II detector, described in detail elsewhere 0], includes a tracking system in 
a solenoidal 1.4 T magnetic field. Electromagnetic and hadronic calorimeters backed by muon detectors surround the 
tracker. The essential detector elements for this analysis are the silicon strip tracking detector (SVX II), the central 
drift chamber (COT) and the central muon system (CMU and CMP). The CMU is a four-layer planar drift chamber 
system outside the CDF magnet coil and calorimeter steel (5 interaction lengths) . The CMP is another muon chamber 
system behind the CMU, shielded by an additional 0.6 m of iron in the flux return yoke. In this analysis, we use only 
information provided by the central sector of the detector, with pseudorapidity | rj |< 0.6. 

Muon candidates are identified by a first-level hardware-based trigger that reconstructs a charged track in four 
axial layers of the COT The trigger then projects the track into the CMU/CMP system and matches the 

projected trajectory to a collection of three or four hits in the CMU muon system within a search window around the 
extrapolated track [q|. The dimuon trigger requires two opposite-sign muon candidates each having px > 1-5 GeV/c. 

The "0(2S) — » A t+ M~ candidates were reconstructed from muon pairs. The ip(2S) events may originate from the 
primary interaction (prompt) or from decays of i?-hadrons (£>-decay). In offline reconstruction, each muon had to 
have at least three hits in the r-<f> strips of SVX II in order to guarantee good vertex information to separate prompt 
and -B-decay candidates. The minimum muon pt is 2 GeV/c. If the CMU candidate has matching hits in the CMP 
chambers [2j, the track px requirement is raised to 3 GeV/c to account for the extra iron traversed. 

The tp{2S) mass and proper time distributions are used in a joint unbinned maximum likelihood fit to extract 
the prompt and B-decay signals in bins of px for ip(2S) candidates. The mass component separates signal from 
background, while the proper time component separates prompt ?/>(2S) events from those produced by B-decays. The 
mass distribution including the radiative tail is described by a combination of a Gaussian plus an asymmetric function 



where mo is a fit parameter for the invariant mass peak, M is the dimuon invariant mass of each event, A is the 
normalization constant, and empirical parameters f3 and n describe the tails of the function. 

The parameters n and f3 of the tail function, and the relative fraction of the Gaussian and CBF are fixed by a fit to 
the entire px range. The Gaussian and the Gaussian part of CBF have the same width a. This width is p^-dependcnt 
due to experimental effects. We use simulation results to describe the relative px dependence in the different bins, 
leaving one width parameter to be fit in the mass probability density function (PDF). The background mass function 



(CBF) given by 
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is linear (P^ ss oc M). 

The proper decay length ct is used to identify prompt and B-decay contributions to the mass signal. Here ct = J^^m > 
where L xy is the transverse decay length projected onto the ip(2S) momentum. The prompt component is described 
by a double Gaussian function centered at zero {Pp 1 )- The long-lived component is an exponential {Pio ng )- Because 
the tp(2S) events from P-decay come from Bh — > ip(2S)X, they do not have a P-hadron lifetime distribution. We 
use the effective lifetime of the S-decay signal as a fit parameter. Because it is defined in the ip(2S) rest frame, it is 
the same in all px bins. Finally, the background in the ct distribution is described by the sum of a prompt double 
Gaussian [P p D plus three exponentials, each convolved with a Gaussian resolution function: one symmetric about 



zero (Pgy m ), one for positive ct only (P£ 

l = f s PT as 

\fsymP sym 



and one for negative ct only (Pf 4 ). The likelihood function is 



[fpPp 1 + (1 — fp)Plong] 
f+Pf + f-P C J + (1 



(l-f s )P^ ss x 



1 



sym 



bkg 
f'+ 



f+-f-)P£] 



(1) 



The population fractions include / s , the ip(2S) signal fraction from the total number of candidates in the fit, f p , 
the fraction of prompt ^>(2S), fsyrm the fraction of symmetric long-lived background, /+, the fraction of positive-ct 
long-lived background, and /_, the fraction of negative-ct long-lived background. The PDFs for the mass fit include 
prnass^ a (Qgjr + Gaussian) distribution for signal, and P b ^ ss . 

The fit parameters are the width of the overall Gaussian in pp ass j the parameters of the double Gaussians for 
the prompt signal and background ct density distributions, the lifetime of the overall long-lived signal, the five data 
fractions from Eq.([T]), and the parameters of the long-lived background functions in each px bin. The fit projections 
in mass and proper time for 5.5 < pr < 6.0 GeV/c are shown in Fig. [1] Projecting the likelihood fits onto the mass 
distribution in each px bin gives x 2 probabilities in the range [0.4 — 100]%. The fitted yields in each px bin are 
summarized in Table |TJ Signal events are classified as prompt or long-lived by the fit parameter f p . For each ip{2S) 
px bin we know the total number of events N. The fit returns the signal fraction f s and its uncertainty cr/ s . The 
signal yield is S = f s ■ N; a 2 s = (aj s ■ N) 2 + N ■ ff. Analogous equations hold for the prompt yield P — f p ■ S and 
the P-decay yield B = (1 — / s ) • S. The correlation between the signal fraction and the prompt fraction is considered 
in the uncertainty of the prompt and P-decay yield. The number of prompt and P-decay events are also listed in 
Table D 
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FIG. 1: The projections on the invariant mass (a) and the proper decay length (b) for 5.5 < pr < 6.0 GeV/c. The fitted curve 
for total signal plus background is overlaid on the data points. The prompt (solid) and long-lived (dashed) proper time signal 
curves are shown separately. 



We have checked the px dependence of all the fit parameters. The variation is smooth and shows no indications of 
rapid changes of the background functions at any px- The prompt fraction decreases approximately linearly in the 
interval 2 < p T < 30 GeV/c. 
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TABLE I: Event yields from the unbinned maximum likelihood fit and prompt event fraction. Transverse momenta are in 
GeV/c. The uncertainties are statistical only. 



The differential cross section is evaluated using the expression 

da(ip(2S)) N(tp(2S)) 



dp T A ■ e reco ■ J Cdt ■ Ap T 



(2) 



Here dg ^^ S - ) - ) is the average cross section for production in the given px bin integrated over rapidity in the 

range | y |< 0.6, N(ip(2S)) is the number of VK2S) events determined by the fit, A is the geometric acceptance combined 
with the CDF dimuon trigger efficiency, e reco is the reconstruction efficiency, J Cdt is the integrated luminosity of the 
data set, and Apx is the width of the px bin. The acceptance and reconstruction efficiency are determined as follows. 

The geometric acceptance is calculated by a Monte Carlo simulation (MC) method, using ^>(2S) — > p + fj,~ decays 
generated uniformly for 1 < px < 40 GeV/c, | y |< 1, and < (j> < 2ir. The ip{2S) decays are handled by evtgen [111 ], 
allowing us to specify the decay polarization as transverse, longitudinal or unpolarized. We generate independent 
MC sets of these three options. Because the tracking proceeds from the large-radius detectors inward, the geometric 
acceptance calculated for the prompt events is insensitive to small displacements of the ip(2S) decay point. Therefore 
we use the same MC samples for calculating the geometric acceptance for both prompt and B-decay events. The 
systematic uncertainty for this assumption is negligible. 

The MC events are passed through the CDF II GEANT-based simulation [3] and the standard CDF reconstruction. 
Events that pass the geometric selection are accepted based on each event's dimuon trigger efficiency, derived from 
CDF data for muon pairs having | y |< 0.6 with each muon having px > 2 GeV/c 13j. Variations with run and 
luminosity arc included in the measurements. The prompt MC sample was analyzed with the likelihood fitter to check 
for pt variations in prompt selection efficiency. None were seen. 

Determining A is sensitive to the "0(2S) polarization parameter a, which defines the muon decay angular distribution 
in the vector meson rest frame: dN/d cos 9 = 1 + a cos 2 9. The polar angle 9 is measured from the vector meson's 
direction in the laboratory frame. 

We have previously measured the ip(2S) polarization in three px bins for prompt events and the average polarization 
in £?-decay events [5|. By symmetry prompt events have a — at px — 0. With 15% probability a x 2 fit shows 
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TABLE II: The differential cross section (pb/GeV/c) times the dimuon branching fraction as a function of pr for | y |< 0.6. 
For the B-decay measurement the symbol Br includes the branching fraction for b quark inclusive decay to V'(2S)X as well as 
the dimuon branching fraction of t/j(2S). 

that the three measured points are consistent with a — 0. We use this as a basis to make the assumption that 
the polarization parameter a is constant over the px range of the data. Averaging the three measured points gives 
an average parameter a = 0.01 ± 0.13, which is used to determine A and its polarization-dependent systematic 
uncertainty. The prompt acceptance A varies from 2% at px = 3 GeV/c to 20% at px — 23 GeV/c. 

For S-decay events we use the same procedure. The polarization dependence is calculated using the measured 
B-decay polarization a e // = 0.36 ± 0.25 ± 0.03 The B-decay acceptance varies from 1.5% at px — 3 GeV/c to 
19% at pr = 23 GeV/c. Since the polarization is different for the prompt and £?-decay events, a weighted average of 
the acceptances in each px bin is used for the inclusive differential cross section. 

The reconstruction efficiency is the product of tracking and muon selection efficiencies measured in CDF data, 
including the tracking efficiencies for the COT (0.996 ± 0.009), SVX II (0.958 ± 0.006) and the dimuon tracking 
and selection efficiency (0.875 ± 0.019). Combining all the factors and adding the uncertainties in quadrature gives 
£reco = 0.805 ± 0.038. 

Because the instantaneous CDF trigger rate might exceed our data handling capacity, the dimuon trigger, like many 
others, is prescaled. The integrated luminosity for the data sample has to be reduced by the luminosity-dependent 
dimuon trigger prescale factor to calculate the cross section. This is done on a run-by-run basis and has negligible 
statistical or systematic uncertainty. The 1.1 fb _1 sample luminosity is reduced to 0.95 fb _1 for this trigger. The 
resulting inclusive cross sections for prompt and i?-decay production are listed in Table HT1 The prompt and J5-decay 
data are plotted versus pr in Fig. [2ja). Data from the Run I CDF measurement [l| are also included in Fig. EJ^b). 

The major systematic uncertainties on these results are due to the systematic uncertainty in the luminosity deter- 
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FIG. 2: (a): Prompt and B-decay production cross section distributions versus pr for these data, (b): The same data with 
the Run I points included. We ignore differences between rapidity and pseudorapidity for this comparison. The B-decay points 
have been scaled down by a factor of ten for clarity of display. 



mination (6%) 1J] and the polarization uncertainty in the acceptance calculation (9% at low px, 2% at high px)- 
Other systematic uncertainties arise from px variations in the trigger (< 3%) and reconstruction efficiencies (4.7%). 
Systematic uncertainties due to the mass shape parametrization, fitting function parametrization, and prompt fraction 
determination are all less than 1%. The data in Fig. [5] have both statistical and systematic uncertainties included. 

The integrated cross sections are calculated by summing the differential cross sections across px bins. The systematic 
uncertainty on the integrated cross section is calculated by assuming that all sources of uncertainty are fully correlated 
among px bins, with the exception of the trigger efficiency uncertainty, which is uncorrelated among px bins. After 
calculating the uncertainty on the integrated cross section from each source, the total uncertainty is calculated by 
summing the individual contributions in quadrature. 

The integrated inclusive differential cross section for px > 2 GeV/c and | y |< 0.6 is measured to be: 



p T (V>(2S)) > 2 GeV/c 



a(pp -> ip(2S)X) ■ Br(?/>(2S) -> pi+ 
= 3.29 ± 0.04(stat.) ± 0.32(syst.) nb. 



For comparison to the Run I measurement, we limit the px range to px > 5 GeV/c. Then the measured integrated 
inclusive cross section is: 



px(ip(2S)) > 5 GeV/c 



a(pp i/j(2S)X) ■ Br(V>(2S) -> ^+p-) 
= 0.69 ± 0.01(stat.) ± Q.06(syst.) nb. 



At 1.8 TeV the integrated inclusive cross section for px > 5 GeV/c and pseudorapidity < 0.6 was 0.57 ± 
0.04^Q Qg nb The increase is (21 ± 19) %, compared to a theoretical prediction of (14 ± 8) % based on changes 
in the parton energy distribution at the higher collision energy 15j. The uncertainty on the experimental ratio is 
dominated by the Run I measurement due to its much lower statistics. 

Prompt ip(2S) production has a harder px spectrum than that for 3/ip production We plot the ratio of these 
two cross sections as a function of vector meson px in Fig. ODJa). The increase in the ratio at larger px reflects 
the slope difference. Even though it neglects any feed-down contributions to J/ip prompt production, the model of 
Ref. Q| predicts the px dependence of this behavior, as shown by the dashed line in Fig. [3](a) . The model prediction 
is normalized to these data in the px bin covering 8-9 GeV/c. This same behavior is seen in the ratio of cross sections 
for production from B-decay, also shown in Fig. [Ha). The ratio of these two ratios is independent of px, as shown in 
Fig-E^b). There is no theoretical motivation for this relation. 



: (b) 
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FIG. 3: (a) The differential cross section ratio of ip(2S) to J/ip as a function of vector meson pr for prompt events (R p ) and 
B-decay events (Rb)- The main error bar on each point is the quadrature sum of the statistical and uncorrelated systematic 
uncertainties for that pr bin. The extensions to the error bars show the correlated uncertainty from the luminosity. The dashed 
line is calculated using the prompt ratio predicted from Ref. 0| normalized to data for the pr bin covering 8-9 GeV/c. (b) The 
ratio of B-decay to prompt ratios Rb/R p as a function of vector meson pr- The fit to a constant gives a % 2 of 13 for 14 degrees 
of freedom. 

In conclusion, we have measured the px dependence of the cross section for ip(2S) production in pp production at 
1.96 TeV. These data have at least an order of magnitude more events than the Run I measurements and show more 
precisely the trends seen in those data. The increase in the inclusive cross section at the higher energy of Run II (1.96 
TeV) compared to Run I (1.8 TeV) agrees with expectations based on the increase in parton energy distribution. 

These data extend the tp(2S) differential cross section measurement up to 30 GcV/c. They are an important input 
for an update of the matrix elements in the NRQCD factorization approach [2]. In the gluon tower model [J], the 
prompt hadroproduction of J/ip, ip (2S), and T states have been calculated. The uncertainties of their calculation 
are rather large but their cross section prediction with adjusted parameters describes the published Tevatron data. 
In addition, their mechanism predicts a longitudinal polarization of J at large transverse momentum which agrees 
qualitatively with the recent Tevatron measurement [5[. We hope that in future calculations with this and other 
models the uncertainties can be reduced in order to make a meaningful comparison to these new cross section data. 
A successful description of both the cross section data and polarization measurements in the perturbative pt region 
would demonstrate a good understanding of the charmonium hadroproduction mechanisms. 
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